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Abstract Ballistic transportation introduces new challenges in the thermodynamic 
properties of a gas of particles. For example, violation of mixing, ergodicity and 
of the fluctuation-dissipation theorem may occur, since all these processes are con- 
nected. In this work, we obtain results for all ranges of diffusion, i.e., both for 
subdiffusion and superdiffusion, where the bath is such that it gives origin to a 
colored noise. In this way we obtain the skewness and the non-Gaussian factor for 
the probability distribution function of the dynamical variable. We put particular 
emphasis on ballistic diffusion, and we demonstrate that in this case, although the 
second law of thermodynamics is preserved, the entropy does not reach a maxi- 
mum and a non-Gaussian behavior occurs. This implies the non-applicability of the 
central limit theorem. 



1 Introduction 



Diffusion is one of the most fundamental mechanisms for the transport of energy, 
mass and information; in physics it is the process by which many systems reach 
uniformity and equilibrium. Due to its importance, it has been the focus of exten- 
sive research in many different disciplines of natural science. The usual manner to 
study the diffusive dynamics is to investigate the mean square displacement of the 
particles, given by (x 2 (i)) oc t a , where (. . .) is an ensemble average. The exponent 
a classifies the type of diffusion: for a = 1, we have normal diffusion; for < a < 1, 
subdiffusion; and a > 1, superdiffusion. The result of this letter applies to all kinds 
of diffusion < a < 2. However, we shall give especial attention to the case a = 2, 
named ballistic diffusion (BD), due to its importance |l|2|3|4j . 

From the point of view of statistical physics, BD has presented some surprising 
features, such as violation of mixing, ergodicity, and of the fluctuation-dissipation 
theorem (FDT) [5 6j. Moreover, ballistic transport is on the boundary between 
stochastic processes described by the generalized Langevin equation (GLE) [7] and 
hydrodynamics. The connection between the violation of ergodicity in the BD and 
elsewhere [8|9| has been discussed recently [10]. Violation of the FDT has also been 
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discussed in systems with activated dynamics [6]. In this context, it is important to 
be sure on what laws, or formalism, one can count on when dealing with BD. Here 
we show that for all diffusive regimes < a < 2, a Gaussian process occurs for the 
probability distribution function (PDF) of the dynamical variable. Furthermore, we 
show that for BD a non-Gaussian process occurs for generic conditions, except if 
the initial PDF is Gaussian. Consequently the central limit theorem (CLT) cannot 
be applied. In fact, this occurs because the system is strongly correlated. We found 
that the entropy grows, but does not reach a maximum, which means that the 
system has energy available for the realization of mechanical work. The proof of 
this is quite general, being valid as long as the noise spectral density is a definite 
even and non- negative function. 



2 Second law of thermodynamics and entropy 



Recent studies |2|4|5|llj have demonstrated that ballistic diffusion possesses peculiar 
properties. Let us consider a Brownian Motion described by a GLE of the form 



dA (t) 
dt 



I T(t-t')A (t') dt' + h(t), 
Jo 



(1) 



where A is a dynamical stochastic variable, T (t) is a kernel, or memory function, 
and h (t) is the noise, which fulfils (h(t)) = and the fluctuation-dissipation theo- 
rem PEG]: 

{h(t)h(t>)) = (A%T(t-t>). 

For an initial distribution of values, A (0), (h (t) A (0)) = 0, it is possible to obtain 
the temporal evolution of the moments of A, 

(A(t)) = (A(0))R(t), (2) 



AHt)) = (A 2 ) p +RHt) (A* (0)) - , (3) 



eq 



2\ 



eq 



where the response function R (t) can be obtained from the inverse Laplace trans- 
form of R (z) = l/[z + r(z)]. For most physical systems, the mixing condition 
R(t —> oo) = occurs. However, this condition fails for ballistic motion (see eq. (J3|)) 
and the time-correlation function of A will be non-null for long times. In other words, 
if the ballistic system is not initially equilibrated, then it will never reach equilib- 
rium and the final result of any measurement will depend on the initial conditions. 
The reason for this is simple, as can be seen from the knowledge that the asymp- 
totic behavior t — > oo can be described by the limit z — > of the Laplace transform 
(final-value theorem). Moreover, it has been shown that if T (z — > 0) oc z v , then 
a = 1 + v [11]. Consequently, for BD, v = 1, T (z — > 0) bz, and 

lim R(t) = lim zR(z) = {I + by 1 ^ 0. (4) 

The major consequence of the violation of mixing is the presence of a residual 
current. Let us suppose that the system starts with an average current Uq such that 



2 



(A (0)} = U . From eq. ©, we obtain (A (*)) = t/o-R (*), where for R (* 00 ) ^ a 
residual current remains. However, the effective current can be very small compared 
to Uq and its value, as any other measurable property for BD, will depend on the 
value of b defined by eq. ©. In other words, the system decays to a metastable 
state and remains in it indefinitely, even in the absence of an external field. Note 
that the "effective friction" 



7 = / T(t)dt = lim bz a - L (5) 
Jo 

is null for all super diffusive motions 1 < a < 2. However, only for ballistic motion the 
major consequences arise. It is important to observe that only for normal diffusion 
with a broadband noise a decoupling for the correlation function 

d -^n-R{t)j\t>)dt> (6) 

corresponding to eq. ([I]) can be done and after a transient time the non-Markovian 
system becomes effectively Markovian. For most of the anomalous diffusions, the 
same can not be done [14] . The major point here is that ballistic diffusion is somehow 
a peculiar kind of motion, being on the borderline between diffusive and hydrody- 
namical behavior. The variance of the stochastic variable A (or simply the mean 
square value of A) can be associated with a temperature by the equipartition the- 
orem. In this way we can write {A 2 (0)) - {A(0)) 2 oc T , where T stands for the 
initial temperature of the system. Equivalently, for a system that reaches equilib- 



With this in mind, eq. ([3]) becomes 



rium, we should have {A / — (A) eq oc T, where T is the reservoir temperature. 



T eff = T + (T -T)[R 2 (t)-l], (7) 

where T e ff is the effective temperature. Notice that some authors [6] have found 
effective temperatures which slowly reach the reservoir temperature after an infinite 
time. In our case, the violation of the mixing condition, R (t — > 00) 7^ 0, implies that 
the system never reaches equilibrium, i.e. T e ff ^ T. 

In order to preserve the second law of thermodynamics, the condition — 1 < R (t — * 00) < 
1 must be satisfied, otherwise heat will flow from the low temperature reservoir to 
the one at high temperature. For that, it is necessary that b > in eq. Now, for 
a system in contact with a heat reservoir (canonical), the memory can be expressed 
as [5E] 

T(i) = J p{uj) cos (wt)dLJ, (8) 

where p (u) is the noise density of states and hni^o F (z) = in BD. From eq. ©, 
we obtain for BD 

6 = lim^= f^>0. (9) 

This simple and general result guarantees that the second law always holds, since 
the spectral density is always non-negative. 
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For a system with specific heat c v (T), one may compute the variation of entropy in 
the process as 

AS = c v (T') (A _ I) dT' > 0. (10) 

Again the second law holds. However, we notice that entropy does not reach a 
maximum since a similar integration from T e jf to the reservoir temperature T will 
produce AS' > 0. Consequently, the final entropy will be 

S = Smax — AS . 

The entropy is not a maximum because the system has not reached equilibrium, 
a common situation in mesoscopic systems |15|16j . which is consistent with the 
violation of mixing and ergodicity, i.e., the system did not forget its initial condition 
completely. 



3 Non-Gaussian behavior 



From the Gibbs entropy, 

S = -k B J P {A) In P {A) dA, (11) 

and Boltzmann's H theorem, the maximum entropy is obtained for a Gaussian 
distribution |17j . Since the entropy is not a maximum, we can guess the distribution 
is not Gaussian. In fact, as will be discussed below, for a BD whose initial PDF is 
not Gaussian it will evolve in the direction of a Gaussian, without reaching it. 

Since eq. (JT]) does not yield the PDF of the dynamical variable A, we need to attack 
the problem using measures of symmetry and non-Gaussian behavior. We can use 
the average values 

(A 3 (t)) = (A 3 (0))R 3 (t) + 3(A(0))(A 2 ) eq 

x[l-R 2 (t)}R(t), (12) 

and 



[A*(t)) = (AH0))R*(t)+3(A 2 ) eq [1-R 2 (t)] 
+2 (A 2 (0))R 2 (t)}, 



(13) 



together with eqs. ([2|) and ([3D, to compute the skewness [18] and the unidimensional 
non-Gaussian indicator [19]. The skewness is a measure of the degree of asymmetry 
of a distribution defined by 

?(*) = { <A 3 (t)) - (A (t)) [3a 2 A (t) + (A (t)) 2 ] } ja\ (t) , 
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where a\{t) = {A 2 (t)) - {A(t)) 2 . Replacing eqs. ©, ©, and {T2]) in q(t), we 
obtain 

,(0)R 3 (t). (14) 
The unidimensional non-Gaussian indicator, defined by 

V {t) = (A*{t))/[3(A* (0> 2 ]-l, 

will determine whether a distribution is Gaussian or not. Replacing eqs. ([3]) and 
([TBI) in the former expression yields 

2 

r?(0)i? 4 (i). (15) 

Equations (|14p and (|15p are important results, since they make explicit the depen- 
dence on the initial conditions. These results are very powerful, being valid for all dif- 
fusive regimes. For example, they state that if the mixing condition R(t — » oo) = 
holds, the final PDF will be symmetric and Gaussian (<r(t) =r)(t) = 0), indepen- 
dent of their initial values. This is true for all diffusions whose exponents a are in 
the range < a < 2. Besides that, if the initial value ?(0) is null then the PDF will 
be symmetric always. In other words, the dynamics preserves symmetry throughout 
the motion. Likewise, if the initial PDF is a Gaussian, then the final one will also 
be. Finally, we have the situation in which the initial distribution is not Gaussian 
(?y(0) 7^ 0) and the mixing condition does not hold. In this case, r]{t — > 00) ^ 0, and 
the final PDF is non-Gaussian; however, it is closer to gaussianity than the initial 
one. This can be seen by analyzing the ratio rj(t) /r)(0) which is always less than 1, as 
long as R 2 (t) < 1. The evolution of the non-Gaussian indicator is in agreement with 
the second law of thermodynamics, since the Gibbs entropy, eq. (jlip . grows as we 
approach a Gaussian. Furthermore, systems that present ballistic superdiffusion are 
strongly correlated, implying a dependence between the stochastic variables A(t), 
due to the action of the memory function during time evolution. With this, we relax 
one of the CLT's conditions, i.e., the condition of independent stochastic variables, 
and we demonstrate that the ballistic transport goes to a non-Gaussian distribution 
depending on the initial conditions. Despite the fact that anomalous diffusion in the 
range < a < 2 is correlated, the correlations are not strong enough, in such a way 
that the CLT prevails and we can impose a relaxation in the independent variables 
condition. 



(JA (0) 
OA it) 



(a 2 (o)) 

{A* (t)) 



4 Numerical applications 



Up to now, we have presented results which are generally valid. Here we shall exhibit 
numerically the mixing violation and the non-Gaussian behavior in BD. First, we 
will begin by defining a spectral density for normal diffusion 

„ 2m < lo < uj d 

P»n{") = { ^ (16) 
, otherwise 
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Fig. 1. Time evolution of the effective temperature (in arbitrary units) for normal 
and ballistic diffusion. We found R(t) numerically for ballistic diffusion (solid curve). 
For BD we use the noise distribution p{uj) = 2jq/tt for 1 < to < 4 and p(u) = 0, 
otherwise. For normal diffusion, we use R(t) = exp(— ■jot). Since BD converges very 
slowly we use 70 = 1 for BD and 70 = 10 -3 for normal diffusion, so we can compare 
both. We choose To = 1.5 and To = 0.5 for the initial temperatures. In both cases 
we used T = 1 for the reservoir temperature. 




Fig. 2. Time evolution of the normalized Non-Gaussian factor rj (t) /rj (0) for normal 
and ballistic diffusion. In both cases, we consider the initial PDF as a Laplace 
distribution with unit mean and unit variance. The functions R(t) are as in fig. ([1]). 

where is the Debye phonon frequency, with 70 being the effective friction, in 
such a way that 7 = lim 2 ^oF(z) = 70, independently of u>d- For ujd/jo S> 1, the 
spectral density approximates a white noise, T(t — t') = 2jo5(t — t') and the process 
becomes Markovian. In this case, R(t) ~ exp(— 7oi). 

Notice that b will be infinite in eq. ([9]) unless we restrict the lower frequency modes. 
In this way, we propose as a noise density a difference between two Ornstein-Zernike 
processes, 

Pbd H = Pw 2 M - Pw x M , ( 17 ) 
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with u>2 > so that, from eq. (jSJ), T (t) = (2jq/tt) [sin (u^i) — sin (u>it)] t . Note 
that in this case lim z _»o T( z ) = 7o — 7o = 0. This null friction explains the existence 
of a residual current, eq. ([2]). A similar spectral density was found numerically in 
the ballistic propagation of spin waves in the disordered quantum Heisenberg spin 
chain [20] . 

Fig. (pQ) displays the effective temperature of a system as a function of time for 
ballistic and normal diffusions. We use eq. ([7]) with T = 1.0 for the reservoir tem- 
perature and with initial conditions To = 1.5 and To = 0.5 for the upper and lower 
curves, respectively. For normal diffusion (dashed curve), we use R (t) = e -70 *, with 
7o = 10~ 3 . For ballistic diffusion we obtain R (t) numerically [14] . using the noise, 
eqs. (fl~6j) and (fT7|) . with jq = 1, u% = 1, and UJ2 = 4. Since BD relaxes very slowly, we 
have to take a "friction" 70 one thousand times larger than that of normal diffusion 
in order to compare both. Notice that, independently of the initial temperature, 
for normal diffusion the system evolves towards thermal equilibrium with the reser- 
voir, converging to its temperature. On the other hand, the temperature of the 
system with ballistic diffusion approaches the reservoir temperature, but does not 
reach it. Instead, it converges to an intermediate value between To and the reservoir 
temperature. 

Fig. ([2]) displays the normalized non-Gaussian factor as a function of time. Here, 
R (t) is as in fig. (pQ). The solid curve displays BD and the dashed curve, normal 
diffusion. Notice that in BD, for large values of t, R(t) oscillates around the value 
predicted by eq. (HJ). For normal diffusion gaussianization occurs, while for BD it 
does not. In both cases, the initial PDF is the Laplace distribution, with (A 2 (0)) = 1 
and (A (0)) = 0. 



5 Conclusion 

We have studied ballistic diffusion for systems with long range memory and stressed 
previous results on the violation of mixing, ergodicity and of the fluctuation-dissipation 
relations. We have shown that the second law of thermodynamics holds; however, 
the temperature does not reach the reservoir temperature and as a consequence the 
entropy does not reach a maximum. Under that condition we show that the velocity 
distribution is not a Gaussian. The study of anomalous diffusion has found many 
applications, from formal mathematical investigation |25|21|22"|2 3 24j to nanode- 
vices [l|2|3|4j . Recently Klumpp and Lipowsky [26] investigated the movement of 
motor particles bound to a cargo particle. They found that the diffusion coefficient 
can be enhanced by two orders of magnitude. Since BD produces less entropy, as 
we demonstrated, and the diffusion coefficient grows with time, it can reach values 
larger than those of normal diffusion. In this way we suggest that an effort should 
be made to design molecular motors based on BD, which should be more efficient 
than those based on normal diffusion. 
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